Spectral scaling laws in MHD turbulence simulations and in the solar wind 
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The question is addressed to what extent incompressible magnetohydrodynamics (MHD) can 
describe random magnetic and velocity fluctuations measured in the solar wind. It is demonstrated 
that distributions of spectral indices for the velocity, magnetic field, and total energy obtained 
from high resolution numerical simulations are qualitatively and quantitatively similar to solar 
wind observations at 1 AU. Both simulations and observations show that in the inertial range the 
magnetic field spectrum Ei, is steeper than the velocity spectrum E v with Et > E v and that the 
residual energy En = Eh — E v decreases nearly following a kj 2 scaling. 
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Introduction. — Plasma motions in astrophysical sys- 
tems are usually magnetized and turbulent. At scales 
larger than characteristic plasma kinetic scales, one-fluid 
magnetohydrodynamics provides a satisfactory frame- 
work for studying such systems [H, 0|. Magnetohydro- 
dynamic turbulence has long been invoked to explain the 
properties of the solar wind, where velocity and magnetic 
field fluctuations are measured in situ over a wide range 
of scales [e.g., UQ- Recent high- resolution numerical 
simulations, however, reported intriguing contradictions 
with the observational data. The Fourier energy spec- 
trum of MHD turbulence obtained from numerical sim- 
ulations appears to have a different scaling compared to 
the scaling inferred from observations. This raises some 
serious questions. Do the numerical simulations correctly 
represent the physics of solar wind fluctuations and, if so, 
then why doesn't solar wind turbulence exhibit the same 
universal scaling found in 3D MHD simulations? 

To formulate the problem, we rewrite the incompress- 
ible MHD equations in terms of the Elsasser variables, 
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that represents a non-dispersive Alfven wave propagat- 
ing along the direction ^va- Nonlinear interactions are 
thus the result of collisions between counter-propagating 
Alfven wave packets. 

Denote by _E ± = (|z ± | 2 )/4 the energies associated with 
the ± waves. Those two quantities are independent inte- 
grals of motion of the ideal MHD system (JTJ. They are re- 
lated to the total energy and cross- helicity, E = E + + E~ 
and H c = E + — E~ , respectively. Cross-helicity pro- 
vides a measure of imbalance between interacting Alfven 
modes; when H c =/= the turbulence is called imbalanced, 
otherwise it is balanced. The solar wind is essentially im- 
balanced, as more Alfven waves propagate away from the 
sun than toward the sun. In a turbulent state, when en- 
ergy is supplied to the system at large scales, both E^ 
cascade toward small scales where they are damped by 
viscosity and resistivity. Theories and numerical simula- 
tions of MHD turbulence address the Fourier spectra of 
the energies E ± (k) in the inertial range of scales, that 
is, scales much smaller than the forcing scales and much 
larger than the dissipation scales. 

According to numerical simulations, statistics of MHD 
turbulence are highly anisotropic with respect to the lo- 
cal mean magnetic field. It has recently been argued 
that the field-perpendicular energy spectra of incom- 
pressible, homogeneous, strong MHD turbulence scale 
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in both balanced and imbalanced 



where the Elsasser variables are defined as z ± = v ± b, 
v is the fluctuating plasma velocity, b is the fluctuating 
magnetic field normalized by y/4npo, v A = B / y/4npa is 
the Alfven velocity corresponding to the uniform mag- 
netic field Bo, P — (p/po + b 2 /2) includes the plasma 
pressure p and the magnetic pressure, po is the con- 
stant mass density, and we neglect driving and dissi- 
pation terms. The guide field Bo can be either im- 
posed by external sources or associated with large scale 
fluctuations that are almost uniform with respect to 
the scales of the inertial range. It follows from these 
equations that for z T (x, t) — 0, an arbitrary function 

z ± (x, t) = F ± (x ± Vytt) is an exact nonlinear solution mismatch in the spectral scalings motivated our interest 



as E ± (kj) c 

cases [fjHlOj, and that this scaling is consistent with an- 
alytic models (llj . This picture seems however to con- 
tradict the observational data of the solar wind, which 
often find the spectrum of magnetic field fluctuations to 
be consistent with the Kolmogorov scaling —5/3 [e.g.,|3j. 
At the same time, recent measurements of velocity fluc- 
tuations in the solar wind reveal an essentially shallower 



spectrum, closer to E v (k) cx k 3 / 2 [e.g., I12h16| . This 
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in the problem. To address the problem we notice that 
there is no strict requirement that the magnetic and ve- 
locity fluctuations be in equipartition with each other. 
Moreover, even though analytic models often appeal to 
the picture of counter-propagating Alfven modes, with 
v = ±b, such Alfven modes are not statistically inde- 
pendent in strong turbulence. This means, in general, 
that (*+■*-) = (v 2 ) - (b 2 ) ^ 0. 

We perform a high-resolution numerical study of both 
balanced and imbalanced MHD turbulence, and con- 
centrate on individual magnetic and velocity spectra. 
We find that in both cases these spectra are gener- 
ally not identical. We observe that the so-called resid- 
ual energy, characterizing the mismatch of the spectra, 
E r (k±) = Eb(k±) — E v (k±), is not universal in that 
its amplitude depends of the driving and the degree of 
imbalance. The scaling of the residual energy is, how- 
ever, close to E r (k) oc kj 2 in both balanced and imbal- 
anced runs. In the balanced case this result was first 
obtained in 0] . While the total energy spectrum is close 
to E(k±) oc fc ± 3 ^ 2 , the presence of residual energy leads 
to steeper magnetic spectrum and a shallower velocity 
spectrum in an inertial interval of limited extent. How- 
ever, since the residual spectrum declines faster than the 
total spectrum, the universality of the turbulence should 
be restored asymptotically at large k±, and it can be 
observed if the inertial interval is large enough. 

For a comparison with the solar wind measurements we 
plot histograms of velocity and magnetic spectral indices 
measured for individual temporal snapshots in numerical 
simulations of MHD turbulence. Comparison of the re- 
sults with analogous histograms obtained from individual 
solar wind measurements reveals good agreement, indi- 
cating that incompressible MHD provides an adequate 
framework for modeling MHD-scale turbulence in the so- 
lar wind. 

Numerical simulations. — The universal properties of 
MHD turbulence are accurately described by neglecting 
the parallel component of the fluctuating fields, associ- 



ated with the pseudo- Alfven mode, e.g., [17H19I ]. By set- 
ting zjj 1 = in equation ((I) we obtain the closed system 
of equations 
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dt TVA - W 



l)- ± + ( 



z T • V_l) z ± = -V ± P 



i/VVS (2) 



in which force and dissipation terms have been added 
to address the case of steadily driven turbulence, and 
we assume that viscosity is equal to resistivity. This 
set of equations is known as the Reduced MHD model 
(RMHD) [mSl]], appropriate for studying MHD turbu- 
lence with a strong guide field, va ^> v rms . Numerical 
simulations of full MHD equations show that the uni- 
versal regime of strong MHD turbulence is reproduced 
well for v a/ v rms > 5 [e.g., 0, H, HH, which is properly 



captured by Reduced MHD system ©. RMHD allows 
one to reduce the number of fields by a factor of two 
and to speed up the numerical integration. We employ 
a fully dealiased Fourier pseudo-spectral method to solve 
equations |(3J) in a rectangular periodic box, with field- 
perpendicular cross section L\ = (2n) 2 and field-parallel 
box size L| = (v a /ty m AL ■> . The choice of a rectangu- 
lar box, as discussed in J1CJ, , allows for the excitation 
of elongated modes at large scales, necessary to avoid a 
long transition region between the forcing scale and the 
beginning of the inertial range. 

The random forcing is applied in Fourier space at 
wave-numbers 1 < k± < 2, fen = 2ir/L». The Fourier 
coefficients inside that range are Gaussian random num- 
bers with amplitude chosen so that the resulting rms 
velocity fluctuations are of order unity. The individ- 
ual random values are refreshed independently for each 
mode on average every t = 0.1 L±/v rms . We define 
the Reynolds number as Re = (L^/2i:)v rms /v. In the 
present simulations, we also introduce correlation be- 
tween v and b, which is achieved by taking as un- 
corrected Gaussian random forces with zero mean and 
variances a 2 . = (|f ± | 2 ). Denoting f v = ^ (f + +f~)> an d 
fb = \ (f + — f~), we obtain that cross-helicity is con- 
trolled by (f„ • f),) = | (er 2 — o 2 _y The results presented 
below have been conducted at numerical resolution of 
1024 2 x 256 points, which corresponds to the Reynolds 
number of Re ~ 5600. In the imbalanced run, the nor- 
malized cross-helicity is H c /E ss 0.6. The cases were run 
for up to 200 large-scale eddy turnover times in order to 
get reliable statistics, see 111. 

The results of numerical simulations are presented in 
Figs. 1 and 2. Two important observations should be 
made here. First, there is a tendency of magnetic energy 
to exceed the kinetic energy at large inertial range scales 
for both balanced and imbalanced cases. The presence 
of nonzero residual energy was noted in previous studies 
[e.g., 0, HH. Comparison of our results with previously 
available numerical data suggests that the level of resid- 
ual energy is not universal, rather, it can be affected by 
the driving and the degree of imbalance. Second, the 
excess of magnetic energy persists in the whole inertial 
interval, however in quite a peculiar fashion. In both 
balanced and imbalanced cases, the residual energy spec- 
trum has a power-law behavior close to E r (k±) cx kj 2 . 
In an inertial interval of limited extent, this leads to 
steepening of the magnetic spectrum and flattening of 
the velocity spectrum, however, the total energy scaling 
stays close to —3/2. Due to the relatively rapid spectral 
decline, the residual energy provides a subdominant con- 
tribution to both kinetic and magnetic energy spectra. 
We therefore propose that the mismatch between veloc- 
ity and magnetic field energies becomes asymptotically 
irrelevant as the inertial range increases, in which case 
the universal scaling —3/2 is restored for both E v (k±) 
and Eb(k±). 
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Spectrum of v and b for balanced MHD 




Residual spectrum E R -E U -E for balanced MHD 



Spectrum of v and b for imbolanced MHD 




Residual spectrum E R =E V — E b for imbalanced MHD 
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FIG. 1: Spectra of kinetic, magnetic, total, and residual en- 
ergies in numerical simulations of balanced and imbalanced 
MHD turbulence. 



Comparison with solar wind data. — Comparison of 
available numerical simulations with solar wind data is 
complicated by the fact that individual solar wind mea- 
surements typically last for a few correlation times (r c ~1 
hour), while in numerical simulations the spectra are av- 
eraged over tens or hundreds of turnover times to ob- 
tain good convergence. To make an appropriate compar- 
ison, we measure the velocity and magnetic field spectra 
for many individual simulation snapshots separated by 
about one eddy turnover time. We then plot distributions 
of the spectral indices obtained in this way for both bal- 
anced and imbalanced cases. The results are presented 
in Fig. (2). 

For comparison, Fig. (j3|) presents analogous histograms 
obtained using solar wind measurements. The data from 
the Advanced Composition Explorer (ACE) spacecraft 
consists of 15,472 different spectra covering the 10 year 
period from 1998 to 2008 [ij]. The data from the Wind 
spacecraft consists of those 120 of 176 spectra stud- 
ied in 14] having the highest normalized cross-helicity 
\a c \ > 0.76, that is, the greatest imbalance. The spectral 
indices are obtained from fits over the range of space- 
craft frame frequencies from 10~ 3 to 3 x 10~ 2 Hz for the 



FIG. 2: Distributions of spectral indices for kinetic, mag- 
netic, and total energies for individual snapshots in numer- 
ical simulations of MHD turbulence. Upper plot: balanced 
turbulence, 80 snapshots; lower plot: imbalanced turbulence, 
196 snapshots. The average spectral indices are indicated by 
arrows. Normal distributions with the mean values and vari- 
ances matching those of the data are also shown. 



Wind data [lj] and from 1.8 x 10~ 4 to 3.9 x 10~ 3 Hz for 
the ACE data It turns out that the scatter of indi- 
vidual indices in numerical simulations closely resembles 
the corresponding scatter in solar wind measurements. 
The excess of magnetic energy over kinetic energy is also 
seen in the solar wind where, in the study [14( for ex- 
ample, the power-law exponent of the residual energy 
takes typical values around 1.75. The solar wind also 
shows a tendency for the spectral indices of velocity and 
magnetic field to be closer together when the normal- 
ized cross-helicity is high than when it is low 0,i. A 
similar tendency is evident in the simulation results in 
Fig. 2. We therefore propose that the mismatch between 
Eb(k±) and E v (k±) observed in the solar wind turbulence 
is neither the manifestation of non-universality of MHD 
turbulence nor does it indicate a breakdown of the ap- 
plicability of incompressible MHD turbulence theory to 
the solar wind. Rather, it is a consequence of significant 
residual energy generated at large scales, in agreement 
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FIG. 3: Histograms of measured spectral indices for the ve- 
locity spectrum (blue triangles), magnetic field spectrum (red 
squares), and total energy spectrum (black circles) in the so- 
lar wind using data from the ACE and Wind spacecraft. The 
average spectral indices are indicated by the arrows. Note the 
different horizontal scales in the two plots. 



with numerical simulations. 

Conclusions. — Velocity spectra, magnetic field spec- 
tra, and total energy spectra in high resolution numer- 
ical simulations of 3D incompressible MHD turbulence 
are shown to be in good agreement with solar wind ob- 
servations at 1 AU where the respective spectral indicies 
of E v , Ef,, and E are approximately centered around 1.4- 
1.5, 1.6-1.7, and 1.5-1.6 [see, e.g., [H 03. It is im- 
portant to note that the large variability found in solar 
wind spectral indices is also observed in temporal snap- 
shots of the numerical simulations. The unique scaling 
laws obtained by averaging simulation spectra over many 
eddy turnover times can also be obtained through a sta- 
tistical analysis of the spectral indices of those snapshots. 
This provides justification for the widely used statistical 
approach to the analysis of spectral indices in the so- 
lar wind, where averaging over many turnover times is 
not practical. Our results indicate that universal inertial 
range dynamics may be present in the solar wind in spite 
of the observed high variability of solar wind measure- 
ments, and that solar wind turbulence spectra are con- 



sistent with the characteristics of incompressible MHD 
turbulence. 
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